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The dynamics of the photoinduced commensurate to incommensurate charge density wave (CDW) 
phase transition in 4ifb-TaSe2 are investigated by femtosecond electron diffraction. In the perturba- 
tive regime the CDW reforms on a 150 ps timescale, which is two orders of magnitude slower than 
in other transition- metal dichalcogenides. We attribute this to a weak coupling between the CDW 
carrying T layers and thus demonstrate the importance of three-dimensionality for the existence of 
CDWs. With increasing optical excitation the phase transition is achieved showing a second order 
character in contrast to the first order behavior in thermal equilibrium. 



Reduced dimensionality seems to be a decisive prop- 
erty governing phenomena like high temperature super- 
conductivity and charge density wave (CDW) formation. 
The latter is typical for quasi one- or two-dimensional 
metals where, in the ground state, the crystal displays 
a static periodic modulation of the conduction electron 
density accompanied by a periodic lattice displacement 
(PLD), both characterized by the wave vector qcDW 
PP. The standard theory that describes the appearance 
of this macroscopic quantum state was formulated by 
Peierls [2]. By considering a one dimensional metal he 
has shown that the divergent static electronic suscep- 
tibility at a wave vector q = 2kp gives rise to an in- 
stability of the electronic system against perturbations 
at this wave vector. This so called Fermi surface nest- 
ing lowers the frequencies of q = 2kp phonons, which 
will eventually evolve into a static lattice displacement. 
From that band gaps at ± kp result, which reduce the 
total electronic energy. If the elastic energy cost to mod- 
ulate atomic positions is lower than the electronic energy 
gain, the CDW state is the preferred ground state. Re- 
cently, this classical picture has been challenged [3H5], 
since (a) the nesting condition derived from the topol- 
ogy of the Fermi surface (2kp) and the observed CDW 
modulation vectors (qcDw) are not generally equal, and 
(b) the diverging susceptibility at q = 2kp is exceed- 
ingly fragile with respect to temperature, scattering or 
imperfect nesting [3 . Contrasting the standard Fermi 
surface nesting scenario the transition from the metallic 
into a CDW state was argued to occur due to strong q- 
dependent electron-phonon coupling [3], particularly in 
transition-metal dichalcogenides [4]. 

Only the concerted interplay of electronic and lattice 
degrees of freedom make the CDW formation possible. 
The two modulations can be individually examined by 
scanning tunneling microscopy [6] , angular resolved pho- 
toemission spectroscopy [5] and electron, x-ray or neu- 
tron diffraction techniques [7 j. In thermal equilibrium, 



the gap in the electronic spectrum and the atomic dis- 
placement amplitude (A) present different projections of 
the same order parameter pp. Adding femtosecond tem- 
poral resolution to the experiment enables investigation 
of their dynamical behavior. Since the electronic sys- 
tem can be perturbed on timescales much faster than the 
characteristic lattice vibration periods these femtosec- 
ond real-time techniques enable studies of the interplay 
between these two components of the order parameter. 
Studies of coherently driven collective modes [8HTD] , pho- 
toinduced phase transitions [TTHT4] as well as manipu- 
lation of the order parameter [TSJ [16] have been per- 
formed with femtosecond resolution in both quasi one- 
and two-dimensional CDW compounds. By tracking the 
photoinduced changes in optical properties or in photo- 
electron emission spectra, information about the tran- 
sient changes in the electronic subsystem of the CDW 
can be obtained, while the lattice dynamics can be de- 
duced only indirectly [10 J. With the recent development 
of time-resolved structural probes, such as ultrafast x- 
ray fT7HT9] and electron diffraction [20H22] techniques, 
structural dynamics in solids are becoming experimen- 
tally accessible [23l [24] , enabling direct insight into the 
coupling strength between electrons and the lattice and 
their dynamical interplay. 

In this Letter we report on systematic excitation den- 
sity studies of structural dynamics in the CDW sys- 
tem 4i^5-TaSe2 in its room temperature phase. Here 
commensurate (C) CDWs exist in the octahedrally co- 
ordinated layers. Upon strong enough excitation the 
phase transition to an incommensurate (IC) high temper- 
ature CDW phase can be achieved on the sub picosecond 
timescale. By means of femtosecond electron diffraction 
we were able to directly probe the structural order pa- 
rameter dynamics following photoexcitation with near- 
infrared femtosecond optical pulses. Resolving the CDW 
suppression dynamics on the sub picosecond timescale 
we demonstrate the importance of specific strongly cou- 



2 




FIG. 1. Crystal structure and charge density waves in 4Hb- 
TaSe2: Panel (a) depicts the alternating stacking of T and 
H layers of the unit cell. The equilibrium diffraction pat- 
tern of the studied sample is shown in panel (b). Here two 
primitive lattice points are indexed. The insert shows the six 
CDW satellites around each Bragg peak. Panel (c): above 
600 K (metallic phase) the Ta ions are hexagonally arranged. 
Below Tc = 410 K the Ta ions of the T layers are commen- 
surately (y/13 x VT5) modulated: twelve ions cluster around 
one central ion to form a commensurate PLD. The conduction 
electron density (shaded) is modulated accordingly. 



pled phonons for the disappearance of the CDW. More- 
over, the CDW reformation dynamics is found to be dra- 
matically slower than in recently studied lT-TaS2 [23] 
and lT-TiSe2 [24] implying the importance of three- 
dimensionality for the existence of CDW order. Finally, 
by studying the non-equilibrium CDW state generated by 
optical perturbation as a function of excitation density 
we demonstrate the change of order of the C-IC phase 
transition: while in thermal equilibrium it is of strongly 
first order the ultrafast optically induced phase transition 
is of second order type. 

4i^5-TaSe2 is a quasi two-dimensional crystal consist- 
ing of covalently bound three atom thick planar layers 
of type Se-Ta-Se which are weakly bound to each other 
along the crystallographic c-axis via van der Waals in- 
teraction, see Fig.[TJa). In the high temperature metallic 
phase 4i^5-TaSe2 belongs to the space group P6s/mmc 
with the hexagonal lattice parameters of a = b = 3.455 A 
and c = 25.15 A [25]. Here, octahedrally coordinated T 
layers alternate with trigonal prismatic H layers, with 
the individual Ta atoms being aligned on top of each 
other along the c-axis. This structure serves as a host 
lattice for CDW induced modulations at lower temper- 
atures: below 600 K an IC-CDW in the T layers is 
formed with qcDW = 0.265 a*. By further lowering 



the temperature through a first order phase transition 
at T c = 410 K a C-CDW with q CD w = 0.277 a* (q CD w 
is tilted by 13.9° from a*) is formed [261127] . The periodic 
lattice displacement of the charge density wave appears 
as (a/13 x a/13)^13.9 modulation, with no modulation 
vector along the c-axis [7]. The modulation, schemati- 
cally shown in Figure [I] (c), is reflected in the electron 
diffraction pattern. Here each of the Bragg peaks of the 
host lattice is surrounded by six weak satellite reflections 
whose intensity is only a few percent of the intensity of 
the Bragg peaks of the host lattice, see Fig.[TJb). In 
both phases, C and IC, no CDW order is present in the 
H layer; a slight modulation therein can be assigned to 
elastic coupling between the two types of layers [7 . Only 
upon further cooling below 75 K, the H layers develop an 
IC-CDW which coexists with the still commensurately 
modulated T layers [7\ [26]. Thus, the CDW in 4H b - 
TaSe2 at room temperature presents a good approxima- 
tion to a true two-dimensional CDW system. 

Single crystals of 4i^5-TaSe2 were grown by chemi- 
cal vapour transport. From the bulk crystal 90 nm thin 
sheets > 150 /im in diameter were cut along the TaSe2 
planes using an ultramicrotome and mounted on TEM 
meshes [28]. Time-resolved electron diffraction exper- 
iments were performed in transmission geometry in a 
collinear excitation scheme applying 200 fs near-infrared 
pump pulses (A = 775 nm) for photoexcitation, and de- 
layed ultrashort electron pulses as probe. We used a 
30 kV femtosecond electron gun with a gold photocath- 
ode, driven by 150 fs laser pulses (A = 258 nm) at a rep- 
etition rate of 1kHz to produce electron bunches con- 
taining ~ 1000 electrons. Their duration at the sample 
position was independently measured by a compact ultra- 
fast streak camera to be 300 fs [29]. The lateral FWHM 
of the electron beam on the sample was 200 /im, while 
its transverse coherence length was measured to be 5 nm. 
Variable time delay between optical pump pulse and 
the electron probe pulse was realized with a computer 
controlled translation stage in the pump laser path. Ac- 
curate spatial overlap of electron and laser pulses on the 
sample was actively controlled by monitoring the position 
of the electron and laser spot and correcting the overlap 
with piezo driven mirrors if necessary. Diffraction pat- 
terns were detected with a pair of chevron stacked micro 
channel plates coupled to a phosphor screen, and imaged 
onto a 16 bit CCD camera. Diffraction images at different 
time delays were taken with 60 seconds exposure (60,000 
shots). Each image with the pump laser pulse incident 
on the sample was followed by an identical unpumped ex- 
posure for reference. The maximum achievable temporal 
resolution was about 400 fs, if < 1000 electrons per bunch 
were used. At large time-delays the data were recorded 
with higher numbers of electrons per pulse to increase the 
signal to noise ratio and decrease the acquisition time. 

We have performed the experiments at room temper- 
ature and analyzed the temporal evolution of photoin- 



duced changes in the intensity of the Bragg peaks and 
their CDW satellites at the individual diffraction orders. 
In all six diffraction orders both satellite and Bragg peaks 
showed identical temporal evolution. Thus, for the sake 
of increasing the signal to noise level, individual CDW 
and Bragg traces were averaged over all available diffrac- 
tion orders [23]. The experiments were performed as a 
function of excitation fluence F, which was varied be- 
tween 1 and 3 mJ/cm 2 . Since the optical beam diameter 
was 500 jam and the optical penetration depth of 110 nm 
[30] is larger than the film thickness nearly homogenous 
excitation is achieved. 

Four diffraction images at different time-delays at 
F = 2.6 mJ/cm 2 are shown in Fig. 2(a). At this ex- 
citation density the intensity of the CDW peaks of the 
C-phase is completely suppressed on a sub-ps timescale, 
corresponding to the photoinduced C-IC phase transi- 
tion. Due to the reduced PLD amplitude in the IC-phase 
the weak CDW peaks of the IC-phase are not resolved. 
Figure [5] (b) and (c) show the evolution of the changes 
in Bragg peak (ieragg) an d CDW satellite (/cdw) inten- 
sities for three different excitation fluences as a function 
of time-delay The suppression of /cdw ranges from 
« 30% at F = 1.7 mJ/cm 2 to a complete suppression 
of the C-CDW order at F = 2.6 mJ/cm 2 . As the CDW 
is partially suppressed, ieragg gains in intensity as the 
crystal is driven towards a more symmetric state, see 
Fig. 2(b). This process is followed by a rapid decrease 
of ^Bragg accompanied by a further reduction in 2cdw> 
which take place on the sub-ps timescale. For the per- 
turbative regime, at F < 2.5 mJ/cm 2 , the recovery of the 
CDW state is observed on the 150 ps timescale, while for 
F > 2.5 mJ/cm 2 the recovery proceeds on a much longer 
timescale. 

In order to analyze the structural dynamics we pro- 
pose the following simple model. Photoexcitation of 
electrons and initial e-e scattering result in pronounced 
changes also in the electron distribution of the Ta 5d 
band. This process gives rise to strong modification of 
the interatomic potential thereby launching highly coop- 
erative (coherent) atomic motion towards the new equi- 
librium positions of the high-T phase via the breath- 
ing amplitude mode of the 13 Ta atom clusters in Fig. 
[l](c). Thus, the PLD amplitude is suppressed as mani- 
fested by the suppression of /cdw and the correspond- 
ing increase in /Bragg- This timescale, T co h? is found to 
be shorter than our experimental time-resolution. Fol- 
lowing the transient increase /Bragg starts to decrease 
on a timescale Ti CO h, during which 2cdw gets further 
suppressed. This timescale is characteristic of electron- 
phonon relaxation processes, whereby the excited q ^ 
phonons give rise to suppression of /Bragg similar to the 
Debye- Waller effect. Similarly, the phonons which are 
strongly coupled to the Ta 5d electrons seem to fur- 
ther destabilize the CDW. Thus, we modeled the rela- 
tive changes of Bragg and CDW scattering intensities 
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FIG. 2. Temporal evolution of the diffraction intensities. 
Panel (a) shows the evolution of the insert in FigJIJb) at sev- 
eral pump probe delays. Panels (b) and (c) represent the 
transient intensity changes in the Bragg and CDW reflections 
at different fluences F, respectively. The dashed lines are fits 
to the data, see text. The vertical solid and dashed lines in- 
dicate the times of the minimum of /cdw and the maximum 
of /Bragg at F = 2.3 mJ/cm 2 , respectively. The inset to panel 
(b) shows the two exponential components which in superpo- 
sition (solid black) give the fit curve for the Bragg transient. 



as Ci(l -exp(-t d /r C oh)) +C 2 (1 -exp(-t d /r iC oh)) J with 
Ci, C 2 < for the case of A/ C dw (*d) an d C\ > 0, C 2 < 
for A/Bragg (^d); the latter has been sketched for clarity 
in inset to Figure |2](b). 

Three modes around 70 cm -1 have previously been re- 
ported to show the characteristic softening of CDW am- 
plitude modes upon warming towards T c [51] . Since the 
ultimate timescale for the (coherent) suppression of the 
C-CDW is given by 1/4 of a period of the amplitude mode 
we chose r co h = 150 fs. Moreover, to fit the entire traces, 
the function has been multiplied by an exponential decay 
to account for the reformation dynamics, and convoluted 
with a Gaussian pulse with FWHM of 400 fs to make up 
for the temporal resolution. From the fit with this sim- 
ple model (dashed lines in Fig. 2) we obtain Ti CO h ~ 500 fs. 
Indeed, these two processes seem to be required to con- 
sistently describe both AIqdw (td) an d A/Bragg (td)- It 
is due to r co h < Ticoh that Icdw reaches its minimum sev- 
eral 100 fs later than the maximum in /Bragg? see solid 
and dashed vertical lines in Figure [2](b) and (c). Inter- 
estingly, this delay is becoming larger with increasing F, 
which can be attributed to different F-dependence of the 
two processes. 

The observation that the complete C-IC phase transi- 
tion is achieved only after a substantial part of the ab- 
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FIG. 3. Fluence dependence of transient PLD displacement 
at td = 1 ps relative to its unperturbed value. The data are 
fit with the BCS function, where temperature was replaced 
by fluence as control parameter. The insert shows the corre- 
sponding changes in icDW- 



sorbed energy is transferred to lattice vibrations via e-ph 
relaxation suggests that the optically induced phase tran- 
sition relies not only on launching the amplitude modes 
via electronically driven change of the interatomic poten- 
tial, but also on exciting specific strongly coupled lattice 
modes. Thus, the photoinduced phase transition seems 
not to be purely electronically driven. Further support 
for the importance of strong e-ph coupling comes from 
the fact that the optically delivered energy required to 
fully drive the C-IC phase transition corresponds well to 
the energy required to heat up the sample to T c = 410 K. 
Using the measured absorbed energy and the specific 
heat of c p = 85.7 J/ (mole K) [32] we indeed obtain a 
temperature increase of ~ 110 K at the critical fluence 
F c = 2.5mJ/cm , see Fig.||| This finding is in agreement 
with previous studies on a quasi-lD CDW in K0.3M0O3 
[T2] and lT-TaS 2 [23 , while in lT-TiSe 2 the absorbed 
optical energy density required to drive the phase tran- 
sition was reported to be lower than the corresponding 
thermal one [24 . 

In the perturbative regime at F < 2.5 mJ/cm 2 the 
CDW reformation is found to proceed with the timescale 
T rec ~ 150 ps, see Fig.[2](c). r rec is about two orders of 
magnitude larger than the observed recovery times in 1T- 
TaS 2 (~ 4ps) [23 and lT-TiSe 2 (« 1 ps) [54]. Thus, the 
CDW in 4i^5-TaSe 2 seems far less cooperative than in 
the two above mentioned dichalcogenides. The most ob- 
vious difference between the polytypes is the interplane 
coupling between the T layers, which is strongly reduced 
in 4i^-TaSe 2 due to the intermediate H layers which are 
not contributing to the CDW. The reduced coupling in 
4i7fc-TaSe 2 is demonstrated also by the reduction of the 
phase transition temperature from T c = 473 K for 1T- 
TaSe 2 to T c = 410 K [27]. In this respect the dramatically 
prolonged CDW recovery in 4i^-TaSe 2 as opposed to the 



other two IT dichalcogenides hints to the importance of 
a certain degree of three-dimensionality for the formation 
of CDW order in these systems. In other words, these re- 
sults suggest that a purely two dimensional CDW order, 
as e.g. in a single T-TaSe 2 layer, might not be possible. 

In the regime above F > 2.5 mJ/cm 2 where the CDW 
is completely suppressed, the CDW recovery proceeds on 
a much longer (> ns) timescale. For this time delays the 
system can be described simply by an elevated tempera- 
ture. Cooling is governed by the heat diffusion, which is 
slow due to the weak thermal coupling of the free stand- 
ing film to the copper mesh. 

To gain additional insights into the nature of the tran- 
sient state, we performed a detailed study of the struc- 
tural order parameter (the displacement amplitude) as 
a function of the excitation density. The displacement 
amplitude relative to the unperturbed state is extracted 
directly from the relative change in /cdw (inset to Fig. 
[3|, taking into account that Iqbw ^ (Q-A) 2 , where Q is 
the CDW wave vector and A is the atomic displacement. 

In thermal equilibrium the C-IC phase transition has 
been identified to be of a strong first order [27]. This 
is consistent with the F-dependence of /cdw in a quasi- 
thermal equilibrium at t^ ^ 1 ns, where a step like change 
between F = 2.3 and 2.6 mJ/cm 2 is observed - see Fig. 
2(c). The functional behavior of the order parameter at 
td = 1 ps upon increasing F (Fig. |3| clearly shows a 
different character. Indeed, the order parameter (|A|) at 
td « 1 ps can be very well fitted by the BCS functional 
form, characteristic for second order phase transitions. 
The dashed line in Fig. [3] is given by a well known ana- 
lytic approximation to the BCS gap equation [33 : here 
instead of temperature T fluence F is used as a control 
parameter, A/A = tanh(1.78(F c /F - 1) 1/2 ), with the 
critical fluence F c = 2.5mJ/cm 2 . This clearly suggests a 
non-thermal character of the photoinduced state. More- 
over, this observation could present an important clue for 
the theoretical understanding of interplay between elec- 
trons and lattice leading to the appearence of CDWs in 
layered dichalcogenides. 

In conclusion, we have studied the structural order pa- 
rameter dynamics in the CDW compound 4i^5-TaSe 2 us- 
ing femtosecond electron diffraction. We show that the 
CDW is suppressed on the sub-ps timescale, launched by 
the electronically driven changes of the interatomic po- 
tential, and further enhanced by strongly coupled lattice 
modes generated via e-ph scattering. In the perturba- 
tive regime the CDW recovery time in 4i^5-TaSe 2 150 
ps) is found to be two orders of magnitude larger than 
previously observed structural recovery times in similar 
transition-metal dichalcogenides. This observation high- 
lights the importance of three-dimensionality for the ex- 
istence of CDWs in layered compounds. Finally, the flu- 
ence dependence of the CDW order parameter showed, 
that the photoinduced ultrafast C-IC phase transition in 
4i^5-TaSe 2 is of second order in contrast to a first order 
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thermally driven (slow) transition. 
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